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We study the = 1/3 quantum Hall state in the presence of random disorder. We calculate the 
topologically invariant Chern number, which is the only quantity known at present to unambiguously 
distinguish between insulating and current carrying states in an interacting system. The mobility 
gap can be determined numerically this way, which is found to agree with experimental value 
semiquantitatively. As the disorder strength increases towards a critical value, both the mobility 
gap and plateau width narrow continuously and ultimately collapse leading to an insulating phase. 
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Two-dimensional electron systems in a perpendicular 
magnetic field have been the focus of both theoretical 
and experimental attention for the past two decades. 
Such systems, if made sufficiently pure and taken to 
low enough temperatureSjexhibit primarily the fractional 
quantum Hall effectdHla (FQHE). While the pure sys- 
tems, especially for the strongest FQHE, are fairly well 
understood, there is essentially no quantitative under- 
standing of the role of random disorder. The impor- 
tance of disorder, however, is underscored by the hall- 
mark FQHE plateaus, which can not occur in a transla- 
tionally invariant system. On the other hand, when the 
disorder strength becomes comparable to the strength of 
interactions between electrons, the FQHE will eventually 
be destroyed. In this paper, we report on finite-size nu- 
merical studies of the effects of random disorder on the 
1/3 FQHE. 

The usual numerically calculated quantities such as the 
energy spectrum, wavefunctions, and the density-density 
correlation functions, while useful in studies of isolated 
impurities 0, 0|, provide little or no understanding of 
transport properties. A more appropriate approach is 
to obtain the quantum Hall conductance by calculating 
topologically invariant Chern integers 0,0,13 in systems 
with periodic boundary conditions (or torus geometry), 
in the presence of random disorder. Physically the Chern 
integers are the boundary condition averaged Hall con- 
ductance of the system in units of e^/h. In the case of 
the integer quantum Hall effect (IQHE), the Chern num- 
bers for the ground states have fixed non-zero integer 
values, while a "metallic" or critical state that separates 
two neighboring IQHE states gives an intrinsically fluc- 
tuating Chern number 8] . (It takes on different integer 
values in response to slight changes of the external pa- 
rameters, such as the disorder configuration). Such a 
quantum phase transition has only been studied for sev- 
eral noninteracting models^, in which only the Chern 
number for single particle states needs to be calculated. 



The situation is very different in the case of FQHE. 
Due to the topological and many-body nature of the 
problem, in the FQHE system at filling factor i/ = 
nh/eB — p/q {B is the magnetic field and n the areal 
density), there exists a manifold of q nearly degenerate 
low energy states on the torus, whose energy differences 
disappear in the thermodynamic limit (Ref . |2| and see be- 
low) . The quantization of the Hall conductance can not 
be tied, in a physically meaningful way, to any particular 
one of the q ground states. For example, in the absence 
of disorder (in which case the degeneracy is exact), when 
an external flux quanta is inserted adiabaticallv|(i llOl 
in a region inaccessible to the electrons, the states within 
a given manifold evolve into each other. In the presence 
of disorder, the Chern number of the individual states 
fluctuates while the sum of the Chern number of all the 
states turns out to be p and robust. As a consequence, 
we may regard the total Chern number p to be shared 
by the q degenerate states, which results in fractionally 
quantized Hall conductance aH = pe'^/qh for the system. 

Based on this picture, we have developed a numeri- 
cal method for studying the topological Chern number 
of the ground state and low energy excited states for an 
interacting system. We will show that the presence of 
disorder leads to several interesting and important re- 
sults for 1/3 FQHE: (i) A weak random disorder lifts 
the degeneracy of the ground state for a finite number of 
electrons. However, the level spacings between the lowest 
three states decrease monotonically with the increase of 
electron number, indicating the recovery of the degener- 
acy in the thermodynamic limit, (ii) The mobility gap, 
which separates the higher energy extended excitations 
from the low energy FQHE states, can be determined 
from the distribution of the Chern number of the many- 
body states. This is the only way to unambiguously dis- 
tinguish between insulating and current carrying states 
in an interacting system, (iii) In general the mobility 
gap, determined this way, will be different from the spec- 
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trum gap (which separates the lowest three states from 
other higher states). It is the mobility gap that should 
be compared with the experimentally obtained activation 
energies, (iv) There exists a critical disorder strength 
W = Wc « 0.2e^/e for Gaussian white noise potential, 
which marks a transition from the FQHE to insulator, 
(v) The physics of the destruction of the FQHE can be 
described as the continuous collapse of the mobility gap; 
the closing of the mobility gap and quantization of the 
Hall conductance occur at the same time. 

We consider a two-dimensional interacting electron 
system in an Li x ^2 square cell with twisted bound- 
ary conditions: T(Lj)<i>(r) = e*^J<i>(r), where T{Lj) is 
the magnetic translation operator and j = 1,2 represents 
the X and y directions, respectively. Calculating the Hall 
conductance au directly from a Kubo formula would re- 
quire a knowledge of all the many body eigenstates for 
a fixed number of electrons. This proves impractical for 
systems with 3 electrons or larger. As first realized by 
Thoulessj^ and co-workers, a topological property of the 
wavefunction[loj . known as the first Chern number, can 
be used to calculate the boundary condition averaged 
an- The importance of Chern numbers, however, goes 
beyond obtaining aH- It appears to be the only quantity 
that distiiiffuishes between insulating and current carry- 
ing states 0, 01 in an interacting system. Because we 
are dealing with many-body wavefunctions, other sim- 
pler numerical methods (such as inverse participation 
ratio and Thouless numbers) used for determining the 
localization of the single-particle wavefunctions have no 
obvious extension here. After a unitary transformation 
= exp[-iY,^J:^{j^Xi + j^y.i)]^k, the boundary-phase 
averaged Hall conductance for the fc-th many-body eigen- 
state is a'^j = C{k)e'^/h, with 

cw.i-/.9,(fe|^)-(^l«}, (1) 

where the closed path integral is along the boundary of a 
unit cell < ^i, ^2 < 27r and summation over j is implied. 
C(fe) is exactly the Berry phase (in units of 27r) accumu- 
lated for such a state when the boundary phase evolves 
along the closed path. To determine the Chern number 
uniquely we separate the boundary phase space into 
approximately 36 — 100 mesh points and get the sum of 
the Berry phases from each mesh. In cases where there 
are near-level-crossings, the integration contour has to be 
chosen reasonably close to these points. This determines 
the size of the mesh, at least locally. For the mesh sizes 
we chose, we found the Chern numbers had converged 
and did not change by further reducing the size of the 
mesh. 

In the presence of a strong magnetic field, one can 
project the Hamiltonian onto the partially-filled, lowest 
Landau level. The projected Hamiltonian in t he p res- 
ence of both Coulomb interaction and disorderllSl can 



be written as: 

H = ^^e-«'/2l/(g)e''i(^--I^^) 

i<j q 

+ T.T.^^'^^'u^^'''''''^ (2) 

i q 

where is the guiding center coordinate of the i-th elec- 
tron, V{q) — /eq is the Coulomb potential, and C/q 
is the impurity potential with the wave vector q. We 
set the magnetic length 1=1 and e^/e£ = 1 for con- 
venience. The Gaussian white noise potential we use 
is generated according to the correlation relation in q- 
space {Uqllq') = [W'^ /A)6q^-q' , which corresponds to 
(C/(r)J7(r')) = W'^6{r — r') in real space, where W is the 
strength of the disorder and A is the area of the system. 
We consider the case oi u = Ne/N^ = 1/3, where 
and Nrj, are the number of electrons and flux quanta. We 
obtain the exact low energy eigenvalues and eigenstates 
using the Lanczos method for systems up to iVe = 8 elec- 
trons in 24 flux quanta, spanning a Hilbert space of size 
A^basis = 735,471. We then calculate the Chern numbers 
using Eq. (1). 

Independent of boundary conditions, we find that, 
for weak disorder, low energy states are separated into 
groups. The lowest group has three closely spaced states. 
This property comes from the three-fold center of mass 
degeneracy of the pure system at = 1/3 and is shared 
by the entire spectrum. For finite sizes, there is a fi- 
nite level spacing between the lowest three states, which 
is much smaller than the energy difference between the 
3rd and the 4th state. The latter is the low energy 
spectrum gap denoted a.s Eg. kt W = 0.06, we have 
Es = 0.04 ± 0.002 (averaged over 100-2000 disorder con- 
figurations), which is size independent for A^e = 5 — 8. It 
remains finite for weak W until W is further increased to 
W « 0.15, at which point Es becomes too small and its 
value in the thermodynamic limit cannot be extrapolated 
from the sizes accessible to our approach. 

On the other hand, the level spacings between the low- 
est three states depend strongly on N^.. They are ex- 
pected to vanish exponentially as the linear dimension of 
the system^]. This can be seen clearly in Fig. 1, where 
we show a semilog plot of the "bandwidth" Ei, vs \/Wl 
for three different disorder strengths. It is apparent that 
E\j drops to zero for large A"e , a direct consequence of the 
topological order in the ground state. This is in contrast 
to the usual effect of disorder in the IQHE, where the 
degeneracy of all the states in a Landau level, not being 
a topological property, will be lifted by the perturbation 
of arbitrarily weak disorder. 

In the presence of weak disorder [W < 0.12), we find 
that the total Chern number carried by the lowest level 
is always 1, for A^e = 3 — 8 and thousands of disorder 
configurations. Note that the vanishing Et and the ex- 
istence of a finite spectrum gap Es, at relatively weak 
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FIG. 1; A semilog plot of the bandwidth Eb of the first level 
that consists of three closely spaced lowest energy states as 
a function of the square root of the electron number N^, at 
three relatively weak disorder strengths W — 0.01, W = 0.06 
and W = 0.12. 
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FIG. 2: The probability distribution P{C) of total Chern 
number C, for the lowest five groups of states. The energy of 
the states increases towards the right. 



W, is a direct manifestation of the = 1/3 FQHE; since 
the lowest three states become degenerate in the ther- 
modynamic limit, the Hall conductance carried by the 
ground state is 1/3 of e'^/h. However, this is not unique 
to the ground state and we find that even the low ly- 
ing excited states above the spectrum gap have similar 
properties. Namely, each three states carry total Chern 
number C = 1. At = 0.01, the lowest two groups have 
C = 1 for all the disorder configurations we have sam- 
pled. In Fig. 2 we plot for the first five groups, P{C), 
for W = 0.06, vs C. P{C) is the probability that the 
total group Chern number is C . We have averaged over 
2000, 1000 and 500 disorder configurations for iVe = 5, 6 



and 7, respectively. As seen in Fig. 2, P(l) = 1 for the 
first group {Ng = 1), and P(C) = for C =^ 1. This 
means all disorder configurations have C — 1, which cor- 
responds to the 1/3 FQHE because each state carries a 
definite average Hall conductance of e^/3/i. For Ng = 2, 
both P(0) and P{2) are nonzero, indicating that a small 
number of configurations have C = or 2. As a result, 
P(l) is reduced to 0.9, 0.91, and 0.92 for Ne = 5, 6, 
and 7. The increase of P(l) with Ne indicates that P(l) 
may recover to 1 at large Ne- For the Ng = 3 case, 
P(l) is significantly reduced to about 0.7; it behaves non- 
monotonically as a function of A^e- This results from the 
coexistence of three different Chern numbers in the ther- 
modynamic limit, which characterizes the delocalization 
of quasi-particle excitations. Namely, these excitations 
carry nonzero Chern numbers which are extended in real 
space. For iV^ > 3, P(l) is further reduced and seems to 
saturate at a value near 0.5. 

We may regard the fluctuation of the Chern number to 
be an indication of the degree of delocalization. In anal- 
ogy to the physics of noninteracting systems we define 
Pext = 1 — P(l) to be the likelihood of the breakdown of 
the quantization of the Hall conductance and thus a mea- 
sure of the delocalization of the charged excitations 0. 
In the FQHE plateau regime, Pext goes to zero as a result 
of the localization or nondissipative nature of the state. 
Here nonzero Pext occurs as we go to higher energy states 
{Ng > 3). The energy that separates these two kinds of 
states is called the mobility edge, where Pext has a large 
increase, which probably becomes a finite jump in the 
thermodynamic limit. For W = 0.06, we find that, from 
Ng — 2 to 3, Pext has the largest increase, which puts 
the Ng — 3 group at the mobility edge. Measuring the 
energy at the mobility edge relative to the energy of the 
lowest level, we get the mobility gap Em for Ag = 4 — 8, 
which is A"e-dependent. For W > 0.06, we determine the 
mobility gap by extrapolating the finite size data to the 
thermodynamic limit. For weaker disorder, the sizes that 
can be treated arc not sufficient to produce a meaningful 
extrapolation to A^ ^ cx3. 

The extrapolated Em vs. W is shown in Fig. In 
the inset we plot Em vs. l/Ag for W = 0.17, which can 
be best fit to Em = O.IO/TV^ + 0.005; thus we obtain 
Em = 0.005 ± 0.003 in the thermodynamic limit. We 
see that at W ~ 0.17 such a gap is strongly reduced, 
consistent with the drop of the spectrum gap for similar 
W, signaling the FQHE is on the verge of being destroyed 
by disorder. 

The energy gap A in the excitation spectrum of the 
correlated many-body ground state can be extracted ex- 
perimentally from the temperature dependence of the 
magnetoresistivity, pxx oc exjp{—A/2kBT), where A/2 
is the activation energy and ka the Boltzmann's con- 
stant Boebinger et al. 0| systematically studied 
the activation energy for i/ = 1/3, 2/3, 4/3, and 5/3 and 
its dependence on sample mobility /i (an indication of dis- 
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FIG. 3: (a) The extrapolated mobility gap Em as a function 
of W. Inset shows at = 0.17 for Ne = 4-8 electrons. 
The dashed line indicates Em can be extrapolated to 0.005 
at — > 0. The dot at IV = is the creation energy for 
a quasiparticle-quasihole pair at infinite separation, extrapo- 
lated from pure systems with up to = 10. (b) Dependence 
of Em on mobility /i. The dashed line is converted from a fit to 
experimental data (taken from Ref . IT3) . Here, we use an em- 
pirical mobility-density relation as well as a mobility-disorder 
relation in the Born approximation (see text for detail). 



order) in GaAs-Ala;Gai_a;As. For a class of high-mobility 
samples, they found that A(/Lt) ~ Co(At)(e^/ei) — F 
consistent with a simple phenomenological model 
that assumes a disorder-broadened excitation energy 
level with half- width F(/i). Figure IHId compares this em- 
pirical formula of A(u) with fitting parameters (Co — 
0.049 and F = 6K [lj|) for high-mobility samples with 
the mobility gap we obtained in our calculation. To do 
this, we assume that both the (zero field) mobility and 
(high field) mobility gap are dominated by short-range 
scatterers (appropriate for high-mobility samples) and, 
in the Born approximation, /i — eft' / {m*^W^) . We also 
use an empirical relation ^ cx n^'^ between /x and electron 
density n (both for zero gate bias) extracted from Fig. 1 
of Ref. Q. Here, we do not include the effects of layer 
thickness and Landau level mixing which, nevertheless, 
exist in experimental samples and are known to reduce 
the gap by as much as a factor of 2 [Hill (for zero or 
very weak disorder). 

As we further increase W, the FQHE becomes unsta- 
ble. This can be discerned by following the evolution of 
the Chern number of the lowest level and an averaged 
over the lowest three states. For example, for A'e = 6, at 
W < 0.14, we have an = e'^/3h; it drops to 0.924eV3/i 
at W = 0.17. At larger W, we find a very strong en- 
hancement of the fluctuation in the Chern number and, 
correspondingly, a rapid reduction of an- Similar re- 
sults are obtained for Ne = 8. As shown in Fig. 4, 
-Ac7h/^W has its largest value near Wc = 0.22 ±0.025 
for all A^e = 5 — 8, which determines the critical disorder 
for the = 1/3 state plateau to insulator transition. 
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FIG. 4: The relative decrease of the Hall conductance of 
the lowest level over the change of disorder strength AW is 
plotted as a function of disorder strength W for A^e = 5, 6, 7 
and 8. 



To illustrate the nature of the phase transition, we 
change the filling to be slightly off 1/3. At weak W, 
the I' — 1/3 FQHE plateau has a finite width due to the 
nonzero mobility gap that survives to fillings slightly be- 
low and above 1/3. But aX W = 0.17, which is slightly 
below the critical Wc, we find the mobility gap _E,„ at 
Ne = Ns/3 + 1 is already reduced to zero. This suggests 
that the plateau for disorder close to Wc is still pinned at 
filling v = 1/3, but with zero width. Therefore, the dis- 
appearance of the FQHE at Wc is caused by the collapse 
of the mobility gap. 
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